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Abstract. Iron-group elements have a very high number of atomic levels and an 
overwhelming number of spectral lines. No NLTE model-atmosphere code can cope 
with these in a classical way. A statistical approach was developed over the last decade 
to decrease the number of levels and lines to a manageable amount. The Iron Opacity 
and Interface (IrOnlc) calculates sampled cross-sections and model-atom files as input 
for model-atmosphere computations. IrOnlc is presently transferred into a parallelized 
code to reduce the calculation time to a reasonable value. It will be accessible by the 
public as a service of the German Astrophysical Virtual Observatory. 



1. Stellar Model-Atmospheres and Statistical Treatment of Iron-Group Lines 

Stellar model-atmosphere codes generally solve the radiation-transfer equation together 
with constraint equations like e.g. hydrostatic and radiative equilibrium. Each of these 
determines a distinct quantity. E.g. the effective temperature T e $ is obtained by solving 
the radiative equilibrium. For NLTE codes, the Saha-Boltzmann equation is replaced 
by the rate equations which leads to a set of coupled, highly non-linear equations. 

A stellar model-atmosphere calculation needs to consider metal-line blanketing to 
create reliable models. Iron-group elements (here Ca-Ni) are very important because 
they have hundreds of thousands of of atomic levels and hundreds of millions of respec- 
tive line transitions. Due to their partly filled 3d and 4s atomic sub-shells, the levels of 
iron-group elements are located in very narrow energy range and, thus, they can be seen 
as a quasi-continuum. 



With classical model atoms (Rau ch & Deetjen 1 120031) for the iron-group ele- 



ments, the numerical capability even of state-of- the-art NLTE codes like the Tubingen 
NLTE Model- Atmosphere Package (TMAfFl. IWerner et al. Il2003h is exceeded by 



far. TMAP has been developed since t he 1980s and is well established (e.g. 
Rauch et al. l2007l : IWassermann et al.ll2010h . It is suitable for hot, compact objects with 



20 000 K 5 T s ff 5 200 000 K, a nd surface gra vities 4 5 log g < 9. It can consider opaci- 
ties from hydrogen to nickel dRauch 112003 ). For the iron-group elements, a statistical 



treatment enables a consideration of an unlimited number of line transitions. 

In order to reduce the quantity of levels and lines for the stellar atmosphere code 
without loosing opacity, IrOnlc has been developed since the late 1990s. The basic 



'http://astro.uni-tuebingen.de/~TMAP 
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Figure 1 . Band structure for Fe vi. The superlevel number is given on the right. 
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Figure 2. Detail of a bound-bound cross-section of Fe vi (band 2 to 6, cf. Fig.jTJ. 



idea of this code was pres e nted b y [Anderson I (1985, 1989) and adopted for TMAP by 
Dreizler & Werner 1 (119921 . Il993h . It combines all atomic levels of an ion into a few 
so-called superlevels. E.g. 241 12 levels of Fe vi are combined into six superlevels 
(Fig-CD- All levels within one superlevel are in LTE relation and contribute to its en- 
ergy and statistical weight. The superlevels are treated by TMAP in NLTE. Transitions 
between these superlevels as well as transitions within the bands are allowed. They are 
calculated with an opacity-sampling technique. In this way, every transition is repre- 
sented with its corr ect strength a nd position (Fig.|2]). Sourc e of the atomic dat a are e.g. 
Kurucz' line lists (iKurucz 1 :2009) and the Opacity Project (Sea ton et al. II 19941) . 



2. The IrOnlc Project 



Presently, the calculation of cross-sections with IrOnlc needs about one to three days. 
As the user cannot start the model-atmosphere calculation until the IrOnlc run is per- 
formed, a new, faster version of this code, using MPI as well as GPU techniques, is 



IrOnlc: How to Consider Hundreds of Millions of Iron-Group Lines 



3 



IrOnlc 

old version new version 



-2 



calculate raw atomic-data values 
create atomic-data files 
create frequency grid 


request via WWW interface 




calculate cross-sections 
elements parallel (manually) 
on different compute cores 


calculate cross-sections 
elements and ions parallel 
on different compute cores 
incl. GPU usage 




create final atomic-data files 






create final atomic-data files 





Figure 3. Scheme of the old (left) and new (right) IrOnlc versions. 



developed. This happens within the framework of a German Astrophysical Virtual Ob- 
servatory (GAVC0) project. 

To speed-up the code, three main approaches will be implemented: 

1. With the old IrOnlc version, the user has to perform three runs to get the atomic 
data file and cross-sections (Fig. 0). The new version will be controlled via a web inter- 
face. The user only has to enter the requested parameters and gets an email notification 
with a wget command to retrieve cross-sections, atomic data, and frequency file. 

2. The new IrOnlc version will use a fixed, fine frequency grid. The cross-sections 
are later interpolated by the model-atmosphere code to the frequency grid of the actual 
calculation. This avoids a complete IrOnlc recalculation in case that the frequency grid 
is changed when e.g. new atoms are taken into account. 

3. The code itself will be transferred into a parallel code. The elements as well as their 
ions will run in parallel with MPI parallelization. The bottleneck of the calculations are 
the millions of theoretical line profiles (approximated by Voigt profiles). Their calcu- 
lation presently consumes nearly 90% of the total calculation time. We tested different 
Voigt function algorithms and got variations up to a factor of ten. Performing these 
calculations on a GPU gives a speed-up of a factor « 40 for the line-profile calculation. 

In summary, a new, faster version of the IrOnlc code to calculate sampled cross- 
sections and model atoms for iron-group elements is written. As this happens within 
a GAVO project, the service will be available to the community. Speed-up will be 
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IrOnlc Request 



Please specify chemical species, line temperature, ionization stages, 
relative abundances (in mass fraction, preset values are solar) as 
well as your e-mail address. 

According to your input, l/Gm/ewH calculate bound-bound and 
bound-free absorption cross-sections based on Kurucz's data 

Once calculated, you will receive a note via your e-mail address for a 
wget of your data. 



Personal Information 



Iron-Group Elements 



Last Name 
First Name 



The model atomsforthe iron-group elements can individually be 
constructed for a single species or for a group of elements combined 
in a so-called generic model atom. 
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Figure 4. The new IrOnlc web interface (presently a test version only). 

achieved through parallelization and restructuring. The new version is estimated to 
be a factor of about 20 faster. This reduces the waiting time until the stellar model- 
atmosphere calculation can start to a value of a few hours instead of days. In addition, a 
web interface (Fig.|4| is created to ensure easy control and accessibility, and to provide 
the data in VO-compliant format. In this way, any model-atmosphere codes can benefit. 
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